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^. . ABSTRACT 

O I Just because the propagator of some field obeys a de Sitter invariant equation 

does not mean it possesses a de Sitter invariant solution. The classic example 
is the propagator of a massless, minimally coupled scalar. We show that 
. the same thing happens for massive scalars with M| < 0, and for massive 

! transverse vectors with My < —2(D — 1)H^, where D is the dimension of 

5^ \ spacetime and H is the Hubble parameter. Although all masses in these 

ranges give infrared divergent mode sums, using dimensional regularization 
(or any other analytic continuation technique) to define the mode sums leads 
to the incorrect conclusion that de Sitter invariant solutions exist except at 
discrete values of the masses. 
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1 Introduction 



Dimesnional regularization |lj is a wonderful tool for perturbative computa- 
tions in quantum field theory but, like all tools, it can be misused. One way 
for this to happen involves the infamous automatic subtraction: dimensional 
regularization registers only logarithmic divergences; it sets power law diver- 
gences to zero. When dimensional regularization is employed to control an 
ultraviolet divergence the automatic subtraction is no problem because the 
right way to deal with ultraviolet divergences is by subtracting them with 
counterterms. In that case the automatic subtraction merely spares one the 
labor of explicitly working out the counterterms required to subtract off any 
power law divergences. However, dimensional regularization can also be used 
to control infrared divergences, and this leads to errors if one fails to recall 
that the technique automatically sets power law divergences to zero. 

The appearance of an infrared divergence in the answer to a quantum 
field theoretic question means that something about the question is unphysi- 
cal. The right thing to do in that case is to revise the question so as to make 
it more physical. The classic example of this is the infrared divergences 
one encounters when computing exclusive scattering amplitudes in quantum 
electrodynamics. Because the photon is massless and all real detectors have 
finite resolution, one can never exclude the possibility that the final state 
contains an extra, very low energy photon. Including arbitrary numbers of 
soft photons in the final state eliminates the infrared divergence [2] . Another 
example is when the vacuum decays, as it does for a massless scalar with 
a cubic interaction. Veneziano showed that even inclusive scattering ampli- 
tudes harbor infrared divergences for this system The problem in this 
case is that the free scalar vacuum centered around = decays, so one 
cannot assume the final state is even stationary, much less centered about 
= 0. Infrared finite results can be obtained by instead releasing the system 
at a finite time in a prepared state centered around = and then following 
its evolution 

A peculiar situation arises in curved backgrounds when parameters of 
the background geometry change the infrared properties of the particle. For 
example, consider a homogeneous, isotropic and spatially flat background, 

ds^ = —dt^ + a?{t)dx-dx . (1) 
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The Hubble parameter H{t) and the deceleration parameter q{t) are, 

The spatial plane wave mode functions for a massless, minimally coupled 
scalar are quite complicated for general q(t) [5] but they take a simple form 
when q(t) is any constant qi, 



The naive mode sum for the scalar propagator between {t, x) and (t', x') is, 

(27r)^-i ^ 
From the small k behavior of the mode functions. 



J k)u*{t', k) + e{t'-t)u*{t, k)u{t', k) \ . (4) 



1<,,<0^,.„- = ^^(^^X-L 1+0(P) . (5) 



0<gi^MM* = — TTT^r^^ 1 X— -n+0(A;")L (6) 



one can easily recognize that the naive mode sum (jl]) possesses infrared 
divergences for all values of qi in the range [H [7] , 

— 1 < gi < (~^^ — ) Infrared Divergences . (7) 

However, most of these infrared divergences are of the power law type which 
are set to zero by dimensional regularization or any other analytic continu- 
ation method. One only encounters logarithmic infrared divergences (either 
in the leading term or one of the k^'^'^+'^'^ corrections) for the cases [7], 

- 1 < gi < =^ qi = ~( ^_2_^2Ar ) for Log Divergence , (8) 

2 

<q ^ qi = { j^_^2N ^ Divergence . (9) 

If one were to incorrectly employ dimensional regularization (or any other 
analytic continuation method) to control the infrared divergence of the naive 



2 



mode sum (jl]) it would seem to give a finite result for the propagator except 
at the discrete qi values ([HHn]) for which there happens to be a logarithmic 
infrared divergence [SI |9j. In fact the mode sum is infrared divergent for all 
values of qi in the range ([7]). The right way of dealing with these infrared 
divergences is not to subtract them but rather to correct whatever unphys- 
ical assumption about the system produced them. In this case the problem 
derives from incorrectly assuming that all modes of the initial state are in co- 
herent Bunch-Davies vacuum, even the ones with super-horizon wavelengths 
which cannot be controlled by a local observer. The system can be made 
infrared finite either by starting with the super-horizon modes in some less 
singular vacuum [TD], or else by working on a spatially compact manifold 
which has no initially super- horizon modes [TI]. Both of these procedures 
augment the naive propagator (jlj) with extra terms which can mediate im- 
portant effects |T2]. 

The purpose of this paper is to point out that a similar situation exists on 
de Sitter background {a{t) = e^*, with H constant) if one considers different 
values of the mass-squared M^. This has important consequences for the 
construction of de Sitter invariant propagators. We show that the formally 
de Sitter invariant mode sums are infrared singular for minimally coupled 
scalars with M| < 0, and for transverse vectors with My < —2{D — 1)H^. 
However, one only encounters logarithmic infrared divergences for the special 
cases, 

Ml = -N{D-1 + N)H'^ and = -{N+2){D-1 + N)H'^ . (10) 

Using dimensional regularization (or any other analytic continuation tech- 
nique) to evaluate the naive mode sums leads to the incorrect conclusion that 
a de Sitter invariant propagator exists, except for the "problematic" cases 
ffTOj) . The correct result is rather that the naive mode sums diverge for any 
scalar with M| < and for any transverse vector with My < —2{D — 1)H'^. 
We make the system infrared finite by working on the compact spatial man- 
ifold T^~^, and we obtain explicit results for the leading infrared corrections 
to the propagator. These corrections break de Sitter invariance, just as has 
long been known occur for the massless, minimally coupled scalar [13] . 

In section 2 we review our conventions for the de Sitter geometry. Section 
3 treats minimally coupled scalars, and section 4 is devoted to transverse 
vectors. Section 5 summarizes and discusses our results. 
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2 The de Sitter Geometry 



We work on the open conformal coordinate submanifold of /^-dimensional 
de Sitter space. A spacetime point x'^ can be decomposed into its temporal 
and spatial components which take values in the ranges, 

— oo < x° < and and — oo < < +00 . (11) 

In these coordinates the invariant element is, 

ds^ = g^udx'^dx" = a^^.Tf^^dx^dx" , (12) 



where rj^y is the Lorentz metric and = —1/Hx^ is the scale factor. The 
parameter H is known as the "Hubble constant" . 

Most of the various propagators between points x'^ and can be ex- 
pressed in terms of the de Sitter length function y{x; z), 



y{x; z) = a^a^H^ 



x — z 



(I 



x^ — z^\—ie 



(13) 



Except for the factor of ie (whose purpose is to enforce Feynman bound- 
ary conditions) the function y{x; z) is closely related to the invariant length 
i(x; z) from x'^ to z^^, 



y{x; z)=4: sin^ (^//£(x; z)J . (14) 

Because y{x; z) is a de Sitter invariant, so too are covariant derivatives of 
it. With the metrics g^u{x) and g^^^z), the first three derivatives of y{x; z) 
furnish a convenient basis of de Sitter invariant bi-tensors [H 



= Ha^(y5l + 2a,HAx,) , (15) 

= Ha,{y5y2a^HAx:), (16) 

d'^y{x; z) 
dx^^dz" 



H^a^a,{y6X + 2a,HAxX-'^aj'^HAx,-2r]^,) . (17) 



Here and subsequently Ax^ = rj^^^x—zy. Acting more covariant derivatives 
just gives back the basis tensors, for example [Ti 



-^^^ = H i2-y)g,.ix) , = ^ (2-?/)^-) • (18) 
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Similarly, the contraction of any pair of the basis tensors produces more basis 
tensors [H], 

= HH4y-y^)=g^^iz)^^, (19) 
- (22) 

3 Scalar s 

The purpose of this section is to demonstrate that the de Sitter invariant 
propagator equation has no de Sitter invariant solution for a scalar with 
< 0, and then to construct the leading de Sitter breaking correction 
terms. We begin by giving the propagator equation and the plane wave 
mode functions for Bunch-Davies vacuum. The associated mode sum can 
be evaluated formally to give a de Sitter invariant hypergeometric function 
which diverges for isolated values of the scalar mass. By studying the infrared 
behavior of the mode sum we show that these isolated values are those for 
which one of the infrared divergences, which are actually present for all M| < 
0, happens to become logarithmic. Except at the isolated values of M|, 
all the infrared divergences are of the power law type which dimensional 
regularization (or any analytic continuation method) incorrectly sets to zero. 
We fix the problem by working on a compact spatial manifold which has no 
initially super-horizon modes, and we derive the leading infrared corrections. 

The propagator of a minimally coupled scalar with mass Ms obeys the 
equation, 

-g{x) [n^ - M|]iA(x; z) = i5^{x-z) . (24) 
The plane wave mode function corresponding to Bunch-Davies vacuum is. 



u{x',k)^^af^' Hl'\-kx') where = ^/ . (25) 
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The Fourier mode sum for the propagator on infinite space is [7], 



(2vr) 



A;)m(2°, k) . (26) 



When this sum exists the result is de Sitter invariant [15 



r(^ + z/)r(^-z/) ^(D-l D-l D y 



(47r: 



r(f) 



2^1 



2 

(4vr)f [^y^ ' '^2 '2 ' 2 '4 



-u: — ; 1 

' 2 ' 4 



ii'^-2r(f-i) r/4N#-i ^ /I 



r(|+^)r(i-^)r(f-i) 



2 ^^2^4^^' (2^) 



r(f)r(i-f 



(27) 



r(| + z/+ra)r(|-z/+ra) .yx „-f +2 



(!)' 



r(i+z.)r(i-z.)„t^o[ r(3-f+n)(n+i)! ^4 

r(^ + z/+n)r(i2^-z/+n)^^y 



r(f +n) n\ 



'■ (29) 



The gamma function r(^ u + n) on the final line of f l29|) diverges for, 



= 1 — ) +N 



= -N{D-1 + N)H^ . (30) 



Its origin can be understood by performing the angular integration in the 
naive mode sum fl26|) and then changing to the dimensionless variable r = 
k/Hy/a^, 



iA'^^(x:z) 



2DTi^H Jo 



dkk^-'^(-kAx) ' j£_3(A;A, 



X 



6{x^-z^)H'^^\-kx^)H';^\-kz''y + ^(2° -x°) (conjugate) , (31) 



X 



(y^r) r) * + ^^(/ (conjugate) | . (32) 

In these and subsequent expressions we define Ax = — That the 
divergence at ( l30l) is infrared can be seen from the small argument expansion 
of the Bessel function and from its relation to the Hankel function, 

Ux) = E \r/f > (33) 
^0 n\r{u+n+l) 

Hi'\x) = ^^^^^^Hiz!^|e-V.(x)-J_.(a;)|. (34) 

The small r behavior of the integrand ( 132|) derives from three factors, the 
first of which is r^~^. The second factor from the Bessel function takes the 
form, 

(-./^^HAxt) ' 'j^f^/E^^HAxr) = —j^J2Ci{ny^ . (35) 
And the final factor from the Hankel functions is, 

H'J^L^r)H<^>L^ry = ^^^^^ £ C.(n).- . (36) 

One does not need the coefficients Ci(n) and C2(n) to see that the small r 
expansion of the integrand takes the form, 

1 oo -p2(,,\n2u oo 

-^12'' A:=0 " ' 1=0 



2V{u)V{2u) 



Y.C,{ny-. (37) 



7rfr(^^)r(z. + l) „=o 

Hence the naive mode sum ( |26l) is infrared divergent for 

- 2 - 2z/ < -1 M| < . (38) 

However, there will only be a logarithmic infrared divergence, either from the 
leading term in (|37j) or from one of the series corrections at n = A^, if one 
has, 

D-2-2v + 2N = -l Ml = -N{D~1 + N)H^ . (39) 



7 



This is precisely the condition f l30p for the formal, de Sitter invariant mode 
sum (12^ to diverge. 

As emphasized in the Introduction, the appeance of an infrared diver- 
gence signals that something is unphysical about the quantity being com- 
puted. The correct response to an infrared divergence is not to subtract it 
off, either explicitly or implicitly with the automatic subtraction of some an- 
alytic regularization technique. One must instead understand the physical 
problem which caused the divergence and then to fix the problem. 

The divergence (138|) occurs because of the way the Bunch-Davies mode 
functions (125|) depend upon k for small k. The unphysical thing about having 
Bunch-Davies vacuum for arbitrarily small k is that no experimentalist can 
causally enforce it (or any other condition) for super-horizon modes. This 
has led to two fixes: 

1. One can continue to work on the spatial manifold R^~^ but assume 
the initial state is released with its super-horizon modes in some less 
singular condition [10]; or 



2. One can work on the compact spatial manifold ^ with its coordinate 
radius chosen such that the initial state has no super-horizon modes 

We will adopt the latter fix. Of course this makes the mode sum discrete but 
the integral approximation should be excellent, and gives a simple expression 
for the propagator which differs from (|26|) only by having an infrared cutoff 
atk = HU 

From the preceding discussion we see that the infrared corrected propa- 
gator iA{x; z) is just (15^ with the lower limit cutoff at r = l/^o^, 

HD-2 .oo JD-3{y/a^HAxT) 
iA(x:z) = pn- / dTT^ ^ ^ 



2^vr— {\^/K^,HAxt)^ 
x|^^(x°-^0)i7«(y^r)/7«(^r)* + ^^(/-x°)(conjugate)|. (40) 



^Making the integral approximation does not alter the renormalization of various = 
scalar models at one loop order [HI [HI [HI [E] , or even at two loops [l9 l (20 1 [2T 1 [22 ] . 
Because the physical graviton polarizations have the same mode functions as scalars with 
M| = |23| . one can also test the integral approximation with the graviton propagator. 
There is no disruption of powerful consistency checks such as the Ward identity at tree 
order and one loop f5T , or the nature of allowed one loop counterterms [251 [ISl [IZj ■ 
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Of course we can express the truncated integral as the full one minus an 
integral over just the infrared, 



dr 



dr 



(41) 

In this case it does not matter if dimensional regularization is used to evaluate 
both iA'^^{x; z) and A^^(x; z) because the errors we make at the lower limits 
will cancel. 

A further simplification is that A^^(x; z) only needs to include the infrared 
singular terms which grow as a^az increases. These terms come entirely from 
the J-„ parts of the Hankel function and they are entirely real, 



(47r)# r(j/+i) 70 



D-2 



"dr T 



J d-3 ( y/a^az HAxt) 



22 



r . (42) 



Before giving the general result for A^^(x; z) it is instructive to work out the 
first two terms in the small r expansion of the integrand. 



r 



D-2 



'1 



y/a^HAxr) 



D-3 
2 



}2u 



ancazH'^Ax'^r'^ 



1 + 



+ 0(t' 



(43) 



Now use the definition (!T3l) of the de Sitter length function to infer. 



^Httr Ha, 



y = ttxttzH^ Ax"^ — 

Hence we have, 

H^-^ 2r(z/)r(2z/) 



axazH'Ax' = {y-2) + {^ + ^) . 

dz dx 

(44) 



A'^^ix-z) 



(47r)f r(^^)r(z/+|) 7o 



X 1 



Z/-2 



2{D-l) ' 2{D-l){v-l)^a, 



■ + - 



r2 + 0(/ 



(45) 
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(47r)#r(i^)r(z.+i 




y-2 , v-{^) la 



2(D-1) 1{D-\){v- 



+---^ (46) 



One can see that the first and second terms of fj46|) respectively cancel the 
= and = 1 divergences in the naive mode sum fl29|) . 
To find the general form of A^^(x; z) we first note (from its expression as 
an integral over fc) that it is annihilated by □ — M|. We next note from ( H6il 
that A^^(x; z) consists of a series of terms, each one of which has the form 
(a^^a^)'^^^^"-'^^ times a series involving powers of {y — 2) and (f^ + f^)- The 
contributions at fixed N must be separately annihilated by □ — , and the 
coefficient of the highest power of {y — 2) at fixed N derives enitrely from 
the A^th order term in the expansion of the Bessel function J ds . These two 

2 

facts imply, 



(4vr)#r(i^)r(z.+ i) 
X E ^ S C-„„„(,-2)"-— , (47) 

where the coefficients CNnm are. 



m\n\iN-n-2m)\ r{^ + N-u) 



ni-u) T(i-u) 



r(^ + A^-2m) r(l-z/+n+2m) r(l-i/+m) ' ^ ^ 

Of course there is no point in extending the sum over N to values N > v — 
-^^) for which the exponent of axttz becomes negative. Those terms rapidly 
approach zero, and they can be dropped without affecting the propagator 
equation because they are separately annihilated by □ — M^. 

We conclude this section by discussing three special cases which occur 
with such frequency as to merit special notation. These are 

Ml = {D-2)H^ ^u = ^ zAb{x; z) = B{y) , (49) 
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M| = ^ i/ 



D-1 
2 

D + 1 



iAA{x;z) = A{y) + SA{ax,az,y), (50) 
iAw{x; z) = + m{a^, a,, y) . (51) 



Although the S-type propagator is de Sitter invariant, its 74-type and H^-type 
cousins have de Sitter breaking parts, 



6 A = kln^Gxaz) , 



SW = J(D-l)Va.-(^)ln(a.a,)(y-2)-(^ + ^)|. 
The constant k is. 



k 



Y{D-l) 



(52) 
(53) 

(54) 



(47r)f r(f) ■ 

The main, de Sitter invariant parts of each propagator consist of a few, 
potentially ultraviolet divergent terms (at y = 0), plus an infinite series, 

H^-"^ { fD N /4xf-i 



B{y) 



+ E 

n=0 



2 

r("-+f ) /2/\n-f+2 r(n+L>-2) /yxn 



(n+1)! M 



(!)' 



r(n+ 



(55) 



(47r; 



4\f-i , r(f+i)/4xf 



-E 



n=l 



r(n + f + l) .y.n-§+2 T{n+D-1) ^y.n 



n-f + 2)(n+l)!M 



(i) 



nr(n+f) M 



(!) 



(56) 



W{y) 



H 



D-2 



'D 



(47r) 



4n#-i 



+ E 



2(f-3)(f-2)Vy 
r(n+f + 2)(f)''-#+^ 



r(f + 2) .4xf-2 
f-2)(f-l)V 

3 ... ... /|/-2^ 

4 

r(n + D)(f)' 



n=2L(^-f + 2)(n-f + l)(n+l)! n(n-l)r(n+f )_ 
And the D-depdendent constants Ai, Wi and W2 are. 



(57) 



A, = 



r(i^-i) 
r(f) ' 



1- 



)+^(^)+^(D-l)+^(l)[, (58) 
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W2 




) -H^) -"'(^+1) (60) 



The infinite series terms of B{y), A{y) and W{y) makes expressions (1551 - 
[57|) seem intimidating. However, note tliat eacli pair of terms in the infinite 
sums cancels for D = 4, so they only need to be retained when multiply- 
ing a potentially divergent quantity. Further, because vanishes more and 
more strongly at coincidence as n increases, only a handful of the smallest n 
terms ever need to be included. This makes loop computations manageable. 
For a massless, minimally coupled scalar with a quartic self-interaction, two 
loop results have been obtained for the expectation value of the stress ten- 
sor [T9], for the scalar self-mass-squared [20] and for the quantum-corrected 
mode functions ^28j. In Yukawa theory it has been used to compute the ex- 
pectation value of the coincident vertex function at two loop order [29], and 
it has been used for a variety of two loop computations in scalar quantum 
electrodynamics [211 [22] . 

The need for de Sitter breaking terms in iAA{x;z) has long been rec- 
ognized [13], and ours reproduces the classic and well known result for the 
coincidence limit of the propagator [SD]. The de Sitter breaking terms also 
show up in the differential equations obeyed by the de Sitter invariant parts 
of the various propagators. 



Whereas the equation for B{y) is homogeneous, the equations for A{y) and 
W{y) both possess inhomgeneous terms which are cancelled by □— M| acting 
on the de Sitter breaking terms SA{ax,az,y) and 6W{ax,az,y). Finally, we 
give some differential relations between the de Sitter invariant parts which 
follow from the series expansions (I55ti57p [31j . 



{Ay-y')B" + Di2-y)B' - {D-2)B 
(%-y2)A" + D(2-y)A' 

{4:y-y'^)W" + D{2-y)W' + DW 




{2-y)A'-k = 25', 
{2-y)W" + ^{D-l)k = 2A" . 



(64) 
(65) 
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4 Vectors 



The purpose of this section is to demonstrate that the de Sitter invariant 
propagator equation possesses no de Sitter invariant solution for a massive 
vector with My < —2{D — 1)H^. We begin by explaining how the full vector 
propagator (longitudinal plus transverse) can be written as the transverse 
vector propagator plus a double gradient of the difference of two known 
scalar propagators. We then derive a formal, de Sitter invariant solution 
for the transverse vector propagator in terms of scalar propagators. One of 
these scalar propagators possesses the infrared divergences we found in the 
previous section. We show how the problem can be corrected, and we derive 
the leading infrared correction. 

The full propagator for a massive vector obeys the equation. 



{x]z) = g^^i5^{x-z) . (66) 



Although this object does appear in certain projection operators, there is 
greater physical interest in its transverse part which obeys. 



D 



[x: z 







D 



[x: z) 



(67) 



An excellent early study of the massive, transverse vector propagator was 
carried out by Allen and Jacobson [32] . A minor error in their work is that 
the source term is not transverse. When this is corrected the propagator 
equation reads [33] . 



-g[x) 



- {D-l)H^ -M, 



A 



[x: z) 



gf,ui5 (x-z) + J-g{x 



d d 
dxf^ dz^ 



i^a{x] z) . (68) 



Although iA^(x; z) contains a de Sitter breaking part, k\n{axaz)i this makes 
no contribution when the propagator is differentiated on both of its argu- 
ments. Therefore, equation (168|) is fully de Sitter invariant. 

Given the transverse vector propagator one can construct the full vector 
propagator by adding the double gradient of a longitudinal part, 



x; z) 



(69) 
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The equation obeyed by iA^{x;z) follows from substituting (1^ in (1^ . 
commuting some derivatives and using fl68p to conclude, 



d d 



iAA(x; z) + 



V 



iA^(x:z) 



0. 



A solution to ([70]) is [31j, 



iA^(x:z) 



Ml 



iA /ifx; z) — iA 



m[x; z 



(70) 



(71) 



^2 - M2 



where iAM{x; z) is the scalar propagator with mass M| — iviy. 

Because we can always construct the full vector propagator from its trans- 
verse part by the procedure just described, we will henceforth concentrate 
on the transverse vector propagator. The most general de Sitter invariant, 
symmetric bi-tensor which obeys the transversality condition f l67p can be ex- 
pressed using an arbitrary function 7(7/) and the basis tensors of section 2 



dS 



(x: z 



{Ay-y')y-iD-l)i2-yh 



+ 



d'^y{x; z) 
dx^^dz" 

(2-y)y-(D-i)7 



dy dy 



■ (72) 



Substituting f l72]) into f pS]) gives the following differential equation for 7(1/) 
away from coincidence (that is, away from x^^ = z^) [33], 

/If 2 n 

(4yV)7" + (^+2)(2-y)7'- [2(Z^-l) + -^]7 = 2(/^-l)5'(y) . (73) 

Recovering the delta function in (168]1 additionally requires that the most 
singular term for ?/ — )• must be [33] . 



i{y) 



H 



D-2 



(74) 



At this point it helps to note that the left hand side of relation f[73]) can 
be expressed as the derivative of a scalar kinetic operator. 

Ml 



{4y-y^)Y + {D + 2){2-y)-f' - [2{D-1) + 

~ IPdy 



7 



{D-2)H^-Ml /[7] , (75) 
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where I[y\ stands for the indefinite integral of 7(7/) with respect to y. Hence 
/[7] obeys the scalar propagator equation, 



{D-2)H'^ - /[7] = 2{D-1)H^B . 



(76) 



This is very like the equation we just solved for the longitudinal part of the 
vector propagator so it should not seem surprising that the unique solution 
for relations fl73ll74l) is, 



7(i/) = 2p-i; 



d r 



M^dy 



E{y) - B{y) 



(77) 



where B{y) is the scalar propagator with mass M| = {D — 2)H^ and E{y) 
is the scalar propagator with mass M| = [D — 2)H'^ + My. To verify fl77|l . 
first note the two leading terms in the series expansions of B{ii) and E{y), 
from expressions (155!) and (129!) with u"^ = (^^j^)^ 



H2 5 



B{y) 
E{y) 



H 



D-2 



(47r)T 
(47r)f 



y 
y 



1 ^fDxfAx%-2 



D 



(78) 



+ ... . (79) 



This establishes that fl77|) has the correct singularity flTl]) near coincidence. 
To check the off-coincidence condition (I75|) . note that B{y) obeys (jHID and 
E{y) obeys, 



{Ay-y^)E" + Di2-y)E' - {D-2)^ 



V 



E = 



Differentiating (16T!) and (IHO!) with respect to y gives. 



(%-y2)5'" + (Z} + 2)(2-y)i?" - 2(D-l)-f 



(%-y2)E'" + {D + 2){2-y)E" - 2(Z}-1)4 



i/2 

Ml 



B' 



E' 



. 



(80) 

(81) 
(82) 



Subtracting (IHT]) from (182!) and multiplying by 2{D — 1);^ shows that our 

ansatz fl77|) indeed obeys equation fl73l) . 

Of course using our ansatz fl77j) in fl72|) only makes sense if the two scalar 
propagators B{y) and Eiy) exist! Recall that they correspond to masses 
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Ml= {D - 2)H'^ and M| = (D - 2)H'^ + , respectively. For D > 2 the 
S-type propagator is safe but the £'-type propagator can have problems if 
My is sufficiently negative. Now consider its formal series expansion f l29|) 

with z/^ = (■^^)^ — ^j^i and in particular, the problematic gamma function, 
V^^^^—u+n). Because it is really the derivative E'{y) which appears in our 
ansatz (1771) . the problem at n = drops out, but there are divergences at 
n = + 1 for all non-negative integers A^. This corresponds to the following 
vector masses, 

D - Is 



+ N + 1 = -{N+2){N+D-1)H^ . (83) 



As emphasized in the previous section, condition (15^ only gives the loga- 
rithmic infrared divergences. For My < —2{D — 1)H^ there will be power 
law infrared divergences which analytic continuation techniques incorrectly 
subtract off. We turn now to fixing this problem. 

Because the infrared divergences we have found for My < —2{D — 1)H'^ 
arise from assuming a de Sitter invariant solution to the propagator equation 
(l68|l . they can only be corrected by abandoning that assumption. How- 
ever, there is no need to discard the formal de Sitter invariant solution 
z[^A^^](x; z). Just like the scalar of the preceding section, we need only 
to add to it a de Sitter breaking, infrared correction. 



x; z] 



dS 



x: z) + 



IR 



[x-z). (84) 



Of course the infrared correction must be symmetric and it must obey the 
transversality condition flB7|) . In analogy with the scalar, we also want it to 
be annihilated by the kinetic operator [□ — [D — 1)H'^ — My]. 

Abandoning de Sitter invariance in [^A[,^](x; 2;) affects the tensor struc- 
ture as well as the scalar coefficient functions. If we preserve spatial ho- 
mogeneity and isotropy (as we did for the scalars of the previous section) 
then the only extra basis tensors we require are the derivatives of u{x\ z) = 
\Yi{axaz) with respect to x^ and z^ 



1 



HaJl , |^ = ifa,5°. (85) 



The most general homogeneous and isotropic tensor with the correct sym- 
metries is. 



^IR 



d'^y dy dy 

{x; z) = F,(a., a., y)^^ + F,{a^, a., y) — — 
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dvi Oy dy On Oti Ou 

The coefficient functions Fi, F2 and F4 must be symmetric under interchange 
of and a^. -F3 does not have this symmetry; when the argument hsts are 
suppressed (to save space) we will indicate the interchange of ax and by a 
superscript T, 

Fj{ax, a^, y) = F^{a^, a^, y) . (87) 
We will also use a prime to denote differentiation with respect to y, 



d 

Fl{ax,a^,y) = —Fi{ax,a^,y) . 



(88) 



The right way to think about the coefficient functions Fi{ax,az,y) is that 
enforcing transversality determines F3 and F4 in terms of Fi and F2. Then 
Fi and F2 are fixed (up to normalization of each independent term) by de- 
manding that [□ - {D-l)H'^ - M^] annihilates [^A^J^]{x; z). 

Covariant derivatives of the new tensor involve some extra identities in 
addition to those of section 2, 



du du 



Dx^'Dx'' " ' dxi'dx'' ' Dz^'Dz 

There are also some new contraction identities, 

du du 



- = -//V(^) 



g^^{x] 



dx^^ dx" 
du dy 



dxf^ dx" 
du d'^y 
dx^^ dx^dzP 



y-2 + 2 

r dy 



.a-,1 



+ 2 



eta 

a. du 



dzP ttx dzP 



du du 

dzP- dz" 
(89) 



(90) 
(91) 

(92) 



The covariant divergence gp''{x)D^i[i,A^^]{x; z) produces one term propor- 
tional to dy/dz'' and another proportional to du/dz'\ The term proportional 
to dy/dzP gives an equation for the coefficient function ^3(0^;, a^, y) in terms 
of Fi{ax, a^, y) and ^2(0^, a^, y), 



dF a dF 



dF dF 

+ {4y-y^)F^ + {D + l){2-y)F2 + (2-y)a,-^ - 2a, . (93) 



da^ 



da^ 
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Because the equation involves only first derivatives with respect to y and a^., 
the general solution can be given using the Method of Characteristics. Once 
^3(0^;, ttzi y) is known, the transversality relation proportional to du/dzP gives 
an equation for the remaining coefficient function Fi{ax-, az,y), 



-i2-y)Fi + {D-l)Fi + ax 



dF, 



I— 

'ax 



'+2^Fi 



dFi a, ^ 

oar At 



+ (iy-y')Fl' + D{2-y)Fi + {2-y)a, 



''X ^x 



dFi 



2a. 



dFi 



.(94) 



dax da^ 

This equation can also be solved by the Method of Characteristics. 

Once one has determined the coefficient functions F^ and F4 from equa- 
tions f l93HM|) . Fi and F2 are obtained, up to normalization of each inde- 
pendent piece, by requiring [^A[J^](a;; z) to solve the homogeneous propgator 
equation, 

Ux - {D~l)H^ - M^] [^^}^] [x- z) = Q. (95) 

Of course this gives a relation for each of the five tensor structures present in 
[^Aj,^](x; z), however, only two of these five relations are independent. The 
one proportional to &^y jdx^dz^ constrains Fi{ax,az,y), 



Ml 



{Ay-y^)F'; + D{2-y)F[ - D + -^ F^ + 2{2-y)F, + 2{2-y)a 



-4a 



dF{ 

' da,.. 



{D-l)a, 



dFi 



dF\ 

da,.. 



d 
da,. 



dFi 



dax 



dax 

2F3 = . (96) 



And -F2(a-X5 o,zi y) is constrained by the term proportional to dyjdx^ dyjdz" ^ 



(Z}+4)(2-y)F^ 
dFl 



-4a, 



da 



2D4 

' {D^\)ax 



Ml 



dF2 



F2-2F[ + 2{2-y)ax 



dF' 



da. 



a,. 



_d_ 

da. 



dFo 



' da. 



dax 

2F^ = . (97) 



The normalization comes from differentiating (!96|) with respect to y and 
then subtracting (1971) to obtain an equation for the combination F[ — F2. By 
making the definition. 



F[{ax, a,, y) - F2(a^., a,, y) = S'{ax, a„ y) 



(98) 



we can identify this relation as the homogeneous equation for an ii^-type 
scalar with mass M^ = {D-2)H^ + M^, 



d 



d 



— (EE]) - (EZD = H-'— Dx - {D-2)H' - Ml 



dy 



dy 



S = 0. 



(99) 



18 



Now write the de Sitter invariant part of the propagator in a form similar to 
the de Sitter breaking part, 



x; z) = ^i{y)^r-7rT, + ^2{y)- 



From relations f l72ll73|) and fl77|) we infer, 

(%-y2)y' + (^ + 2)(2-y)7'-2p-l)7 



(100) 



^( - ^2 = - 
Hence we infer 



4(L)-l)if2 



E'{y) 

2m 



(101) 



F[{a^, a^, y) - ^2(0^, a^, y) = 

1 H^-^ r(z/)r(2z/ 



J-^A'^'ix-z) , 
2H^dy ^ ' ^ ' 



(102) 



2m (4^)f r(^)r(z.+|) ^„ ^^-(^)-iv 

xE (- + -)" E (iV+l-ri-2m)C^+i_(y-2)^— ^-,(103) 

n=0 '^z '^a; m=0 

where A^^(x; 2;) was defined in equation (H7|) . the coefficients Cjynm are given 
in fllg]) . and the index z/ obeys z/^ = (^^^)^ — 

It is instructive to give the first (A^ = 0) term in the series expansions of 
the four coefficient functions Fi{ax,az,y), 



iF,)o 



7 



(y-2) 



(z/+^)(z/-l) 



D+l ' 



D+l ' 



D+1- 



+ — ) + ( — 



D-3^ 



(z/-l)^a2 a^;^ ^a^ 



a.xflz 



(104) 
(105) 
(106) 



D+l \ n,, D+3 



V- 



(z/-l) ^a^ a 



+4- 



D-3\ 



V- 



0+1 

2 



.(107) 
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Here the index v obeys 
7 



'■^2^)^ — -j^ and the constant T is, 



1 E 



D-2 



TT 2 



r 



T{v)T{2v) 



D+l- 



(108) 



2 ' 2^ 

The general series expansion for Fi{ax., az, y) has the form, 



N=0 



V — 



-N 

n=0 



a. 



m=0 



Af+l-n-2m 



(109) 



As with the de Sitter breaking corrections for the scalar, there is no point in 
extending the series for N > u — The other coefficient functions have 
the same first line as (11091) so we give only their subsequent forms. 



r iV— n i 



F2 



m=0 



^ n 



2m 



(110) 



n=0 



a. 



m=0 L 



p3a f^x 0.2 
^ Nnm\ "T 



"T ^ Nnni\ 

dz 



(1/-2) 



N-n-2m 



N+1 



r iV+2~n i 



E(7^+?)" E 



Nnm 



(Z/-2) 



N+2-n-2m 



;iii) 



(112) 



n=0 



m=0 



We close this section by giving the transverse vector propagator for the 
important special case of. My = —2{D — 1)H'^. For that mass the associated 
scalar has M| = —DH^, corresponding to the W^-type propagator considered 
at the end of section 3. For this reason we give the vector propagator a 
subscript W, and we decompose it into a de Sitter invariant part and a de 
Sitter breaking part. 



At 



ix: z 



[x:z) + 



x; z) 



;ii3) 



The de Sitter invariant part is, 

{Ay-y^){W-Br + {D-l){2-y){W-By 



(x; z) 



A{D-l)m 



d'^y 
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+ 



{2-y){W-By' + {D-l){W-By '\ dy dy 
A{D-1)H^ \dx^d? 



V 



.(114) 



And the de Sitter breaking part is, 




-ln(a.a,)(y-2) + Z}(^ + ^) 




V 




ln(a:i.a^) 



9y dy ttx du dy 





.(115) 



5 Discussion 

We have shown that infrared divergences preclude de Sitter invariant so- 
lutions for the propagators of either a minimally coupled scalar with mass 
< 0, or for a transverse vector with mass My < —2{D — 1)H'^. (If one 
includes the longitudinal part of the vector propagator then infrared diver- 
gences occur for My < 0.) However, in most cases these infrared divergences 
are of the power law type which is automatically subtracted by any regu- 
larization which is based upon analytic continuation. (We stress that these 
considerations apply as well to the standard technique of continuation from 
Euclidean de Sitter space.) Only the special values of M| and My given 
in equation f lTOj) result in logarithmic infrared divergences which show up 
in analjd;ic regularization techniques. Thus one might reach the incorrect 
conclusion that de Sitter invariant propagators exist for all scalar and vector 
masses, except for a few "singular" cases. 

That conclusion is wrong because infrared divergences should not be 
renormalized the way one treats an ultraviolet divergence. The appearnace of 
an infrared divergence signals that an unphysical assumption has been made, 
and the right response is to identify the problematic assumption and modify 
it. In our case the unphysical assumption is that the universe can have been 
prepared in coherent Bunch-Davies vacuum for arbitrarily long wavelength 
modes. There is no causal process by which this can be accomplished in the 
de Sitter geometry. When one assumes either that the initially super-horizon 
modes are in some less singular state [10], or else that the spatial manifold 
is compact [11], the resulting propagators become infrared finite, but not de 
Sitter invariant. 
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In each case the true propagator can be written as the naive, de Sitter 
invariant result (defined by dimensional regularization) plus a de Sitter break- 
ing infrared correction which is real and obeys the homogeneous propagator 
equation. For the scalar propagator we have, 



iA(x:z) = iA'^^(x:z) + A 



x: z 



(116) 



where iA'^^{x; z) is expression (1291) and A^^ (x; z) is given by relations 
The analogous expression for the transverse vector is. 



[x: z 



dS 



,x: z) + 



IR 



;ii7) 



where the de Sitter invariant part is defined by expressions (!72|) and (!77|) . and 
the de Sitter breaking terms are given by equations and (1108111121) . The 
full vector propagator, including the logitudinal part, is given by equations 
dM]) and dni). 

It might be wondered what physical sense it makes to consider the prop- 
agators of particles with tachyonic masses. First, there is the mathematical 
point that they don't possess de Sitter invariant propagators, despite what 
one might conclude by erroneously defining these propagators with some 
analytic continuation technique. Second, there is the important issue of fol- 
lowing the time dependent vacuum decay which must occur when symmetry 
breaking takes place during a phase of infiation. In this respect the infrared 
correction terms may be quite important, as they sometimes are for the anal- 
ogous case of FRW geometries with constant deceleration [12] . 

Another application for our propagators is the projection operators for 
higher spin propagators, in which case there are no physical particles with 
tachyonic masses to worry about. For example, consider the graviton h/j^u in 
exact de Bonder gauge. 



pa 







;ii8) 



Just as the source term for the transverse vector propagator equation (168|) 
must be consistent with Lorentz gauge (167|) . so too the source term of the 
graviton propagator equation must be consistent with (11181) . The resulting 
projection operator turns out to involve the full vector propagator. 



V 

fiu ' pa 



[x]x') = gp{pga)uiS^{x-x') - -^--^g^,u9paiS^{x-x') 



1 

+ 2 



D^D'ptlA^]{x; x') + D^D; z[,Ap](x; x') 
D,D'tlA,]{x;x') + DM.tlAp]{x;x') 



.(119) 
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One can easily check that the de Bonder gauge condition, 



- \D'Jf"'] % LPpJ (x; x') = . (120) 



2 

requires the vector to have mass My = —2{D — 1)H'^ 



g{x) [n^ + {D-1)H^\ I [^A,J (x; x') = g^.id'' {x-x') . (121) 

This is not only tachyonic, it actually corresponds to the first of the special 
cases (!83l) for which the transverse part harbors a logarithmic infrared di- 
vergence, so the problem would show up even in an analytic regularization 
technique. 

This is all highly relevant to the debate concerning the de Sitter invari- 
ance of free gravitons |3l]. It has long been obvious to cosmologists that 
free gravitons cannot be de Sitter invariant because they share the same 
mode functions as massless, minimally coupled scalars [23]. Indeed, the ten- 
sor contribution to the primordial anisotropics of the cosmic ray microwave 
background derives from precisely the same infrared singular dependence of 
these mode functions [35]. On the other hand, some relativists insist that free 
gravitons must be de Sitter invariant because de Sitter invariant solutions ex- 
ist for the propagator when a gauge fixing term is added to the action [36] . 
In previous work we have shown that adding these gauge fixing terms is not 
valid, owing to the linearization instability [211 EZ]- Imposing an exact gauge 
condition such as flllSp should still be all right, but we have just seen that it 
leads to an inevitable breaking of de Sitter invariance through the projection 
operator. (This same breaking can be seen as well in noninvariant gauges 
by adding the appropriate compensating gauge transformations [38].) Anto- 
niadis and Mottola long ago discovered a similar problem in another gauge 
[39]. Before our current work one might have dismissed these examples as 
"spurious IR divergences for the Feynman propagator in the sense that the 
IR divergences are absent for other values of gauge parameters" [IQ]. They 
now appear as just those cases for which infrared divergences, that are al- 
ways present and which always break de Sitter invariance, just happen to go 
from being of the power law type to logarithmic, and hence become visible to 
analytic continuation techniques. And the correct procedure is not to ignore 
them or subtract them but rather to remove the erroneous assumption of de 
Sitter invariance. 
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